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摘 要 : 利用 积分 恒等式 和 插值 后 处 理 技术 ， 本 文 在 各 向 异性 网 格 上 对 Sobolev 型 方程 的 Carey 非 协调 有 
限 元 解 进行 高 精度 算法 分 析 。 首 先 ， 根 据 Carey 元 的 特性 ， 即 其 有 限 元 解 的 线性 插值 和 线性 元 解 
相同 ， 我 们 构造 插值 后 处 理 算 子 ， 得 到 了 有 限 元 解 的 超 盘 近 性 质 和 整体 超收 敛 及 后 验 误 差 估计 。 
接着 ,根据 误差 渐 近 展开 式 ， 运 用 外 推 方法 ， 进 一 步 得 到 了 具有 四 阶 精度 的 近似 解 。 
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有 限 元 的 插值 后 处 理 、 外 推 、 校 正 等 方法 是 提高 有 限 元 解 精度 的 有 效 方法 ， 它 们 已 被 广泛 
地 应 用 于 计算 的 理论 与 实践 中 tt" 引 。 其 中 文献 [1-3] 对 协调 元 的 超收 敛 、 外 推 等 进行 了 研究 。 文 
献 [ 稀 研究 了 非 协 调 Wilson 元 的 超收 敛 性， 对 于 一 类 二 阶 椭圆 边 值 问题 ， 得 到 了 一 类 特殊 网 格 
下 ( 即 要 求 所 有 单元 的 水 平 边 和 竖 直 边 分 别 相等 )， 和 矩形 Wilson 元 在 顶点 和 四 边 中 点 的 超收 敛 
性 质 。 而 文献 {四 基 于 Wilson 元 的 渐 近 误差 展开 ， 对 于 二 阶 椭圆 边 值 问题 得 到 了 整体 超收 敛 和 
局 部 超收 和 敛 以 及 校正 格式 。Carey 元 是 一 个 著名 的 四 自由 度 非 协调 三 角形 元 ， 关 于 它 的 收敛 性 
分 析 已 有 许多 研究 结果 (参见 文献 [6-11])。 其 中 文献 [8] 给 出 了 Carey 元 在 正则 网 格 下 关于 二 阶 
椭圆 问题 的 整体 超收 敛 、 外 推 和 校正 格式 。 基 于 文献 (12] 中 的 各 向 异性 插值 定理 ， 文 献 [9| 与 
文献 [10] 研究 了 Carey 元 对 二 阶 椭圆 问 题 的 收敛 性 和 超收 敛 性 ， 并 用 数值 算 例 验证 了 理论 分 析 
的 正确 性 。 文 献 [11] 则 把 它 元 应 用 到 各 向 异性 网 格 上 位 移 障 碍 下 二 阶 变 分 不 等 式 问 题 ， 并 得 到 
了 最 优 误差 估计 。 

考虑 下 列 Sobolev 型 方程 


-Au —-Au=f 在 9x(0,7] 中， 
u=0 在 aN x (0,T] E, (1) 
u(z,y,0) = g(x,y) 在 Q, 


其 中 Q C RR 是 有 界 矩 形 区 域 ， 其 边界 99 分别 平 行 于 zx 一 轴 和 y 一 轴 ，f e IL2(Q), g 是 一 个 光 
滑 函数 。 其 等 价 的 变 分 问题 为 : 求 u€ Hg (Q) 满足 下列 方程 


| (Vur, Vv) + (Vu, Vv) = (fjv, Yve HLO), 


u(0) =g. 


(2) 
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JR (1) 广泛 地 存在 于 许多 实际 物理 过 程 中 。 许 多 作者 用 有 限 差分 和 有 限 元 法 对 此 问题 进 
行 了 讨论 出 二， 而 文献 | 引 给 出 了 其 双 线 性 元 的 整体 外 推 ， 文 献 [1 各 给 出 了 其 Wilson 元 解 的 高 
精度 分 析 ， 得 到 了 Oh) 阶 的 整体 超收 敛 和 后 验 误 差 估 计 ， 以 及 O(h9) 阶 的 外 推 结 果 ， 但 上 述 
所 有 分 析 均 基于 网 格 的 正则 性 假设 。 本 文 将 利用 文献 [2,15] 中 提出 的 高 精度 方法 ， 并 基于 文 
RR [10,12] 的 研究 导出 各 向 异性 Carey 元 解 O(h2) 阶 的 整体 超收 敛 和 后 验 误差 估计 。 最 后 根据 误 
差 渐 近 展 开 式 得 到 了 比 文献 14] 高 一 阶 ( 即 O(h4) 阶 ) 的 整体 外 推 结 果 。 


2 Sobolev 型 方程 的 Carey 元 逼近 及 整体 超收 敛 分 析 


假定 Jn 为 Q 的 一 族 直角 三 角形 前 分 ， 对 任意 的 天 € 几 的 两 条 边 分 别 平行 于 z- 轴 和 ?一 轴 ， 
并 且 所 有 水 平 边 和 竖 直 边 分 别 相等 ， 但 这 里 不 要 求 满足 正则 条 件 和 拟 一 致 假设 ( 记 为 GATM 三 
AJEPI UA), V^ c mo) 为 相应 的 线性 有 限 元 空间 ， 并 求 密 e Vp. E 


| (Và, Vv) + (Và^, Vv) = (fov), Vve Vf, ài 


a^ (0) a Ing. 
这 里 以 及 以 后 出 现 的 1, 为 线性 插值 算 子 。 
变 分 问题 (2) 对 应 的 Carey 70383129: ok u^ e VE， 使 得 
as (u*,v) -- a(u^,v) = (fv) VveVg, 
u^(0) — Ing. 


其 中 
an(w,v) = > I Vw.Vvdzdy, Ing € Vj. 
KEJa K 
iki Xd (4| HAREE Wilson 元 的 点 态 超收 敛 时 曾 用 到 GATM 网 格 ， 文 献 [2] 中 用 到 的 
等 腰 直 角 三 角形 网 格 显然 是 GATM 网 格 的 一 种 特例 。 但 文献 [2] 和 文献 身 都 要 求 剖 分 满足 正则 
性 假设 。 
由 文献 [2] 和 文献 [14 可 得 下 面 引 理 。 
引 理 1 3€ J4, 9 GATM 三 角形 网 格 剖 分 ，w fl a^ 分 别 为 Carey 元 和 线性 元 的 解 ， 则 
Ipu” = ^, Inu? =}, 


BI Carey 元 解 的 线性 插值 和 线性 元 解 相 同 。 
定理 1 Ru, 记分 别 是 (2) 和 线性 元 空间 的 有 限 元 解 ， 则 


t 
le^ - nul < C1?[ f (lulis lula) ?as] (5) 
0 


证 明 ”由 方程 (2) (3), XER [10] 的 引 理 1 及 Cauchy 不 等 式 可 知 ， 对 任意 的 ve V7 


(VU — Inus), Vv) + (V(ü^ — Inu), Vv) 


(Vu, — Inu), Vv) + (V(u — Inu), Vv) 


^ 


< Ch? (julla + lulls)lols < Ch (uelis + llis)" + [vfi (6) 
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在 (6) 式 中 取 v = (G — nu) I 


p 2 ld,4 2 2. | 2 
|a 一 ul 十 2di ü^ - Inul < Cht (lulls + llulla) + |a 一 Inu p (7) 


BR 
dja 
Fl -Thuli < CR (lulls + lulls)? 


XpEXPSRE. JHEXSIS"(0)— Ihu(0)， 可 证 本 定理 。 
为 了 取得 整体 超收 敛 和 后 验 误差 估计 ， 设 天 & 四 是 由 相 邻 四 个 及 中 的 单元 合并 构成 的 一 
个 大 单元 ， 我 们 构造 插值 后 处 理 算 子 [19] 
IL, «cao > RE). 


定理 2， 设 IT2 为 上 述 构造 的 插值 后 处 理 算 子 ， 那 么 在 引 理 1 的 假设 下 ， 有 
2 eR 2 i 2 2 2 
ITE, - s], «en CE f eti + tas] tas. 
i-a = [e TE], «005. 
证 明 注意 到 


2 h EEG. ED 2 7 2 7 
u^ —u- [] u^—uz][] u^ — T] u [| u—u. 
2h guo 2h Ph an PY t loin 


由 引 理 1 和 定理 1 
| IL, Inu — IL, Iall, 


t i 
= [IL a^ - n], «eta -iru er [ f ua rutas] , 


再 由 插值 基本 误差 估计 
Tm =o, = TE, » < niai. 
可 得 | l 
HIE, ^ ul, «ee (Cf ate al)as] + hes} 
同 理 可 得 l 
Iv — uli = fe - TL, l|, +00. 
3 ME 


引 理 2 UJ, X GATM 三 角形 网 格 剖 分 ，w e H*(Q), Hve V^, Wl 


f VwVvdzdy = E / (hzuzzy — hyuzy,) v; dzdy 
Q 12 Jo 


hs 
-u J (hyuzyy — hztzzy)vydzdy + O(h®)lulslvh. 
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证 明 ”同文 献 [16|， 设 + = 二 KKU Ki, ?7-KUK;. BM cuc 


h 
n = -2 f (hzUzrzy 一 hyuzyy) vzdrdy T O(h^)|uls,- lvli, 


为 此 考察 函数 


P "E 1 » MEOS 
B(ü,9) = I (@ — In) bedédn + F (Geen — tem) Pedédn. 


根据 Sobolev RA EM, IER E H*?(7) M9 c Vh, RIE 
IB(9, 9)| < CIl elo,- 


直接 计算 可 得 
B(àü,v) = (, Vüc PT), V Ve € PaP). 


根据 Bramble-Hilbert 8| 38 
IBG, P) € Clélaslblos, V 9 € P). 


由 可 逆 映 射 Fk : K UK KU Ky 和 积分 变换 可 得 (8) 式 。 
同 理 可 证 
h 


n E f (hyuzyy — hzuzzy)vydzdy + O(h*)|uls,. |v]1,. 


12 


最 后 由 (8) 和 上 式 ， 对 所 有 单元 求 和 即 可 证 明 引 理 。 
引 理 3 在 五!- 模 意义 下 ， 成 立 


Inu? = Ipu + RO -O(h*), Or e b. 
证 明 0 — U^ Inu, H32, XHEX vc ob, 


(V6, Vv) + (V0, Vv) = (Vu — Inu), Vv) + (V(u — Inu), Vv) 


其 中 
h 
万 = T [hs (uina 本 Uszy) 一 hy (utzyy * Usyy)]: 
hz 
h= 12 [hy (uzyy + Uzyy) — ha(utzzy + uzzy)]. 
考虑 下 面 辅助 问题 


一 Au — Aw = ns (fis + fag) 在 Qx (0, T] 中 ， 
w=0 在 aQ x (0,T] L, 
w(x,y,0)=0 fr Q P. 


-f frosdady — | fonydrdy Oi) (lulls + lulls) lvla: 
Q 


(10) 
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Xu u E€ H5(Q)， 则 有 广 € H1(Q)， 那 么 由 微分 方程 解 的 正则 性 理论 有 
lloc lla + ollis < Clwells + fills). 
上 述 辅 助 问题 (10) 的 等 价 变 分 问题 为 : 求 we HEA) 
| (Vus, Vv) + (Vo, Vv) = — fo fivsdrdy — fo favydrdy, Vv € HEO), a 
w(0) = 0. 
设 we Hi(0) &o^ e Vh 分 别 为 辅助 变 分 问题 (11) 的 准确 解 与 有 限 元 解 ， 则 有 
(V(&, — 33^), vv) + (V(6 — h20), Vv) = O(h*) (lluis + lulls} loh: 
在 上 式 中 取 v —6,— oh, ur 
所 8 一 Pet 上 < Ch? (lusis + lulls). 


两 端 积 分 并 注意 到 0(0) 一 h2@*(0) = 0， 即 可 证 明 本 引 理 。 
利用 和 定理 1 相同 的 方法 可 证 下 面 引 理 。 
引 理 4 Bw, o^ 4 3 (11) 的 准确 解 和 线性 有 限 元 解 ， 则 


t 1 
o^ - e|, < cé [ f. (alls + ol)’as] 
0 


为 了 获得 整体 渐 近 展开 ， 我 们 构造 一 个 Jan 上 的 Lagrange 四 次 插值 后 处 理 算 子 TIw, 。 由 插 
值 基 本 误差 估计 和 插值 的 性 质 以 及 文 [16] 的 方法 可 知 ， 在 GATM 三 角形 网 格 剖 分 上 ， 对 任意 
的 we 五 5(O)， 下 面 三 式 成 立 。 


4 4 4 zx 
= 2 h, 
| NICE, as | | Lel, < Cl|wli, Yw ev; 


IIT, e-e], < lols 07324. 
定理 3 设 .J 为 GATM 三 角形 剖 分 Jan 中 点 加 密 得 到 的 网 格 ，w*(Y = L3) AE 
F Jyn (y = 1, 3) W Carey 元 解 
W ||u — v^ |;  C(u)h*. 
证 明 dd CT ID 5 MOX SIS 及 引 理 4 知 


4 4 D to É 
4h uv -u- IL, (n — Inu — KO") +h? IL, o" + 


= hw e (TE, 9^ -w) + On) 


Il 


hw + h2 IL," sd) v ID, o - uw) + (14) 


h2w + O(h*), 
即 
IL. u^ — u 4. hw 4- Olh’). 
于 是 经 一 次 Richardson 外 推 可 得 : u^ = w 十 O(h4)， 定 理 得 证 。 
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Higher Accuracy Analysis for the Anisotropic Carey Element 
Solution to Sobolev Type Equation 


SHI Dong-yang!, HAO Xiao-bin? 
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Abstract: By using the integral identities and the interpolation postprocessing technique, the higher 
accuracy approximation of the anisotropic nonconforming Carey element for solving the Sobolev type 
equations is investigated. Firstly, the interpolation operator is constructed, the superclose. global 
superconvergence and posteriori error estimate are obtained with the help of the distinct property of 
Carey elements, i.e., the linear interpolation of the solution for Carey elements is equal to the solution 
for linear triangular element. Secondly, by virtue of the extrapolation method, the accuracy of the 
related approximate solution with fourth order is derived through the asymptotic error expansion. 
Keywords: Sobolev type equations; Carey element; higher accuracy 


